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Abstract 

From the epoch of recombination (z ~ 10 3 ) till today, the typical den- 
sity contrasts have grown by a factor of about 10 6 in a Friedmann universe 
with fl = 1. However, during the same epoch the typical gravitational po- 
tential has grown only by a factor of order unity. We present theoretical 
arguments explaining the origin of this approximate constancy of gravita- 
tional potential. This fact can be exploited to provide a new, powerful, 
approximation scheme to study the formation of nonlinear structures in 
the universe. The essential idea of this method is to evolve the initial dis- 
tribution of particles using a gravitational potential frozen in time. We 
carry out this scheme for several standard models including the CDM and 
HDM and show that the results match quite well with those obtained by 
exact numerical simulations. We compute different statistical measures of 
clustering and compare them for the description of nonlinear evolution. 

This approximation also provides valuable insight into understanding 
various features of nonlinear evolution; for example, it provides a simple 
explanation as to why pancakes remain thin during the evolution even in the 
absence of any artificial, adhesion-like, damping terms. We also compare 
this approximation with other schemes like Zeldovich approximation and 
frozen-flow. Our procedure has a far greater range of validity than the 
Zeldovich approximation since it can handle motion across (and inside) a 
caustic properly. Unlike in frozen-flow, actual shell-crossing does occur in 
the frozen-potential approximation; hence it provides a far more accurate 
description of the velocity field compared to frozen flow approximation. 
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1. Introduction 



It is believed that structures like galaxies and clusters of galaxies 
formed out of small inhomogeneities via gravitational instability. In this 
paper we discuss a new approximation scheme called Frozen Potential Ap- 
proximation (hereafter FPA) ( Bagla and Padmanabhan 1993a ; Brainerd 
et. al. 1993) for studying the nonlinear growth of the inhomogeneities. In 
section 2 we introduce and compare the various approximation schemes and 
in section 3 we discuss the results obtained from numerical simulations. In 
section 4 we present some analytical results which can be obtained using 
FPA, and the limits of applicability of FPA are studied in section 5. 



2. Quasilinear Approximations 

In studying the process of structure formation, linear perturbation the- 
ory can be used when the density contrast 6 is small but one has to rely 
on full N-body simulations to grasp the nonlinear evolution. Our under- 
standing of processes in the nonlinear regime is very limited due to the 
absence of analytical results for the nonlinear phase. However, there exist 
some approximate schemes for evolving the trajectories of particles which 
are valid in the quasi-linear phase. These approximations are very useful 
in understanding the physical processes in the nonlinear regime. 

The exact trajectory of a particle in a Friedmann universe is described 
by the equations 

_ + 2-- = -^V* V^ = 4*<W* (1) 

where a(t) is the expansion factor, x is the comoving coordinate, and p b is 
the background matter density. We shall set Q = 1 and consider pertur- 
bations only in the matter dominated era, though the results can be easily 
generalised for other cases. These equations can be recast using scale factor 
as the time variable, and by introducing the variable if) — (2/3H 2 )<f>. In 
terms of a and i/>, (1) and (2) become 

d 2 x 3 cfx 3 „ , „ , 6 

— + = V^; V 2 ^ = - (2) 

da z 2a da 2a a 

As we shall see, various approximations follow very naturally from 

these equations. 
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It is apparent from (2) that the velocity of particles is not constant but 
decays even when Vt/j — due to the damping term (3/2a)(dx/ da) in (2). 
In Zeldovich approximation (hereafter ZA) ( Zeldovich 1970 ; Shandarin 
and Zeldovich 1989), it is assumed that the velocity of a particle remains 
constant along its trajectory and has the value it had at the starting point. 
This is equivalant to taking 

d 2 x dx „ , . 

so that the trajectory in ZA is 

X = x *n-«WL ( 3 ) 

By assuming a uniform velocity, in contrast to the decaying velocity 
in the absence of a potential, some amount of accleration has been taken 
into account. This approximation works well till the trajectories of parti- 
cles cross but leads to unphysical results after crossing. To see how this 
happens, consider the motion of particles near a potential well. All the 
velocity vectors in the neighbourhood of a minima of the potential point 
towards it and all the particles move towards the minima. However, af- 
ter the particles cross the minima, they do not turn around and fall back 
into the minima, as they should in a self gravitating system of particles. 
Instead, the particles continue travelling with their initial velocity leading 
to dispersal of structures. Hence this approximation works well only if the 
average displacement of particles is smaller than the scale over which Vt/j 
changes sign. 

Zeldovich ansatz can be used to show that first structures to form in 
a generic gravitational collapse are planar surfaces of high density, the so 
called pancakes. The unrealistic part of dynamics in ZA is that particles 
after collapsing and forming a pancake continue to move with the same 
velocity leading to thickening of pancakes. 

The Frozen flow approximation (FFA) ( Matarrese et. al. 1992 ) over- 
comes the thickening of pancakes by assuming that the motion of particles 
is well approximated by potential flow. In FFA the inertia of particles is 
assumed to play a negligible role in comparison with the expansion of the 
Universe or a rapidly ( spatially) varying potential. Therefore the first term 
in (2) can be neglected and the motion may be described by 
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d> !X! d>"3£. , , , , , , 

^ = 0; -(x, a) = -W(x,0 (4) 

The gravitational potential at the initial time generates a velocity field 
which is used to move the particles. Particles reach the minima of potential 
only asymptotically and the pancakes do not thicken, since (strictly speak- 
ing) no pancake has really formed. The dynamics in this approximation 
is unrealistic near pancakes as it does not predict oscillations of particles 
about the pancakes, a generic feature of motion in such a case. 

It is possible to generate a better approximation (Bagla and Padman- 
abhan, 1993a; Brainerd et al., 1993) by the following procedure : we know 
that the gravitational potential is a constant in the linear regime for an 
Q — 1 Universe. It can be argued ( based on models like spherical top 
hat) that the potential does not change significantly even in the nonlinear 
regime. This approximate constancy of potential can be used to evolve den- 
sity perturbations to mildly non-linear density contrasts with much more 
realistic dynamics. Frozen Potential Approximation (FPA) assumes that 
the potential remains constant even in the nonlinear regime. The trajecto- 
ries are then evolved using (2) but with a fixed potential: 

cf 2 x 3 cfx 3 „ , . 

By using the full nonlinear equation, a realistic dynamical evolution of tra- 
jectories is obtained. In FPA, near a minima, the particles move towards 
the minima of the potential and then oscillate about it with decaying am- 
plitude due to the expansion of the Universe. In this approximation, shell 
crossing does occur and the pancakes are formed; the particles oscillate 
about the pancake with a small amplitude and move along the pancakes 
towards the minima of the potential. 

The trajectories of a pair of particles in a potential well is plotted for 
the three approximations ( ZA, FFA, FPA ) in Fig.l. The potential is 
taken to be ij) = — cos(x) and the trajectory of particles is plotted on the 
x — a plane. In Zeldovich approximation, the particles move with uniform 
velocity and even after reaching the minima of potential, continue travelling 
with the same velocity, erasing the structure in the process; the trajectory 
in this case is therefore a straight line ( dot - dashed ). In FFA, the 
particles approach the minima only asymptotically ( dashed line ), giving 
the impression of a stable structure even though no structure has really 
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formed in finite time. On the other hand, FPA gives trajectories oscillating 
about the minima with a decaying amplitude ( unbroken line ). In exact 
N-body simulations the damping is more prominent and the pancakes are 
thinner than the ones obtained by FPA, because FPA does not take into 
account the deepening of potential wells as particles cluster. 



D 
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Fig.l : Trajectories of particles in a one dimensional potential well in various 
approximations are shown in the x — a plane. The dot-dashed line corresponds to 
ZA, the dashed line to FFA and the continuous line to FPA. Notice that in FPA, 
the particles cross the potential minima and oscillate about it. 
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3. Numerical Simulations 



We tested three approximations using numerical simulations for mod- 
els such as the CDM and HDM and compared the results with exact N-body 
simulations. In the simulations, we start with an initial potential ij) which 
is a realisation of a gaussian random field with a power spectrum corre- 
sponding to the model under consideration PAk) — P s (k)/k i . The power 
spectra were normalised by using the COBE quadrapole value. The poten- 
tial was obtained on the lattice sites of a cubic (128 X 128 X 128) grid. The 
particles were moved by the ansatzs described above. 

Various approximations can be compared visually in Fig. 2 and Fig. 3 
which give slices [ with dimensions (96 x 96 x 15) (/i -1 MPc) 3 ] of the Uni- 
verse in these schemes at different epochs for CDM and HDM respectivly. 
It is obvious that at early stages {z 1) of structure formation, all approx- 
imations give similar results, while at later stages (z ~ 0) the differences are 
quite prominent. Zeldovich approximation is the worst of the three; after 
shell crossing, the particles continue moving with the same initial velocity 
leading to thickening of the pancakes. In FFA the particles are deposited in 
the pancakes and then they move slowly towards the minima; the pancakes 
are very thin and there is very little dumpiness even for CDM spectrum. 
Pancakes do not thicken in FPA and, at late stages, clumpy structures 
dominate. 

In Fig. 4 we have shown the density contrast er(r) defined as a 2 = £ 

with 

?(r) = 4 f r ti*)d* 3 = ^P- 

for three redshifts for standard CDM using FPA. For comparison we have 
also plotted the N-body result at z = 0. We have obtained the N-body 
result by using the universal scaling relation between the linear density 
contrast and the actual nonlinear density contrast (Hamilton et.al. 1991 ; 
Nityananda and Padmanabhan 1993). In the linear regime the evolution 
is self similar, i.e. the density contrast grows in proportion with the scale 
factor. When high density contrast is reached, FPA tends to underestimate 
the density contrast. 
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The FPA can be quite useful in studying the evolution of voids in 
the Universe. To compare the sizes of voids with theoretical results it is 
necessary to choose a suitable definition of void. One way of estimating the 
sizes of voids is to compute the rms displaceent of particles in a simulation 
(Shandarin, 1992) which can be related to the average diameter of voids by 
D vold — 2kd rms where k is a constant of order unity which depends on the 
shape of voids. For spherical voids k = 1.3. The thick line in Fig. 5 shows 
the evolution of d rms for CDM and HDM. For the standard CDM model, 
we found that d rma ~ 18/iij~ MPc at z = 0. 
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Fig. 5 '. Evolution of d rms . Here d rms is plotted as a function of the scale factor, 
the continuous line is for CDM and the dashed line corresponds to HDM model. 
Dot-dashed lines have been drawn for the power law solutions, a and a 1 / 1-8 for 
early and late times respectively, as discussed in the text. 



4. FPA : Analytical approach 

The equation (5) can be solved analytically in some special cases (Bagla 
and Padmanabhan, 1993b). Consider, for example, the motion of particles 
near a local density minima. The density contrast very close to the extrema 
is a constant and can be represented by 8 — —A with < A < 1 a constant. 
Since for constant density, V?/> — — rA/3, (5) becomes 

d 2 r 3 dr _ A 

da? + 2a"da ~ ~a ' ' 



7 



Bagla and Padmanahhan 



This has the solution 

r = ^£ =s i n h(v / 2Aa). 

Since the density contrast is valid for small r and hence for small a ; in this 
limit, r oc a, as expected from linear theory. This suggests that the radius 
of voids increases linearly with the scale factor for a <C 1 • 

The density profile away from the centre of the void can be approxi- 
mated as 

where n is a positive constant. This leads to the equation 
d 2 r 3 dr 3 L n -, „ 

1 = r 1_n (81 

da 2 2a da 2a (3 - n) J 

This nonlinear equation has one simple solution : 

1 l/n 



3n 2 



La x l n . (9) 



(3 - n)(n + 2) 
which holds for n < 3. 

For a gaussian random field, the profile of density around a density 
extrema is a power law with the same index as the correlation function( 
Bardeen et. al. 1986 (BBKS)). For our Universe, the index of the correla- 
tion function is about 1.8, therefore we expect the voids to grow as a 1 / 1 ' 8 
or (1 + .a) -1 / 1 ' 8 in the late stages. In the linear limit, we expect the void 
radius to grow as a. While calculating d rms from simulations, we average 
over regions with different index n and hence expect d rms to grow at a rate 
intermediate between a and a 1 / 1 ' 8 . At late times merging of voids will slow 
down the growth rate even further. 

These features are shown in fig 5. The two dot- dash lines indicate 
growth proportional to a and a 1 / 1,8 and have intercepts chosen to match 
with d rms for CDM. It is clear that these lines match with d rms at early and 
late times respectively. Since the degree of nonlinearity is lower in HDM 
compared to CDM (when normalised to COBE), the d rms oc a for HDM in 
this range. 

Similar analysis can be carried out for generic motion near a density 
peaks and it can be shown that very close to the peak, the solutions are 
oscillatory. 
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5. FPA : Accuracy and Limitations 



The results above show that FPA is a good approximation for a CDM 
spectrum and other similar spectra. The FPA is based on the assumption 
that the potential ip changes very little during the evolution of density 
perturbations. To understand the agreement with N-body simulations and 
to study the limitations of FPA, we should ask the question: how much 
does the gravitational potential change in reality ? 
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Fig. 6 : Variation in Potential with scale factor. We have plotted the variance of 
•potential £ as a funtion of the scale factor for r = 1,5, 10, 50 and 100/i _1 MPc. 
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At large scales, at which a 2 < 1 and linear theory is valid, the gravita- 
tional potential is a constant to a high degree of accuracy. Similarly at very 
small scales - dominated by fully virialised structures - the gravitational 
potential does not change significantly. The evolution of gravitational po- 
tential in the quasi-linear regime depends on the power spectrum. It is 
possible to obtain a semianalytic form for this evolution by using the form 
of the scaled relative pair velocity h = (v patr / Hr) where v pair {x) is the 
average relative pair velocity at the scale x and r — ax. From the N-body 
simulations, it has been observed that h(x,t) — /i[<t 2 (:k, t)] is a univer- 
sal function of a 2 (Hamilton et al., 1991; Nityananda and Padmanabhan, 
1993). We can use this universal relation to obtain a relation between the 
linear density contrast and the nonlinear density contrast. We find that the 
relation between true density contrast a 1 and the linear density contrast 
a\ is well approximated by: 

((T 2 L (l,a) (for cr 2 < 1 ) 
<r 2 (x,a) oc 1 k7 2 (/,a)] 3 (for 3 £ a 2 £ 50) (10) 
I a)] (for 50 <> a 2 ) 

where/ = x(l-\-a 2 ) J / 3 (Bagla and Padmanabhan, 1993b). If we approximate 
<r^(/, a) oc a 2 l~( n+3 ^ where n is a local index, this gives us for (<r/a) : 

a(x,a) I a (»+ 4 ) 10 < a 2 < 50 

OC < („ + 2 ) (11) 

a [a-Wi 50 < g 2 

In the standard CDM model n ~ —2 in the nonlinear domain (cr 2 <; 50) 
and n ~ (—1 to —2) in the quasilinear regime. From (11) we see that (c/a) 
remains constant for n = — 1 and n — —2 \n the quasilinear and nonlinear 
regimes respectivly. Hence in CDM like spectra, there is a conspiracy of in- 
dices ensuring that there is very little change in the gravitational potential. 
The above analysis, of course, is approximate since the relation between i\> 
and (a /a) is nonlocal and the index n varies with scale. The true variation 
of tp can be obtained by integrating the Poisson's equation numerically. 
The Fig 6 shows ip at four different scales in standard CDM model, as a 
function of redshift. 

The FPA, however, is not accurate at small scales, for a different rea- 
son. The location of pancake like structures in our ansatz is influenced 
significantly by the initial potential. This implies that the large scale 
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structures hardly move in our ansatz. In reality if the streaming veloc- 
ity at a (large) length scale L is v(L) then structures at this scale would 
have, on the average moved a distance of H^ 1 v(L). In CDM models at 
scales L ~ 50/1" 1 Mpc, v ~ 200& ms : hence these structures could have 
moved about H^v ~ 2h~ 1 Mpc by now. Thus at scales smaller than about 
\h~ l Mpc we would expect the approximation to show some inaccuracy. 

These aspects can be clearly seen in fig. 7a, b which compares FPA 
with exact N-body simulations in 2D. The pancakes in N-body simulations 
are thinner than those in FPA (due to deepening of potential wells) and 
are also slightly shifted with respect to those in FPA. 

6. Discussion and Conclusions 

The FPA is a powerful approximation which can be used to obtain 
approximate results, numerically as well analytically. It compares well with 
the N-body simulations in a statistical sense. By construction, the minima 
of potential do not move in this approximation and this makes it inaccurate 
at small scales. A very strong point in favour of this approximation is that 
it contains full information about velocities and therefore it can be used to 
study various models in the redshift space. The FPA can also be used to 
study the motion of baryons in dark matter potential wells. 
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Figure Captions (For figures not included here) 
Fig. 2 : Evolution of density perturbations in CDM for various approx- 
imations. These frames, from left to right, correspond to a = 0.25, 0.5 
and 1.0 ( present epoch) respectively. The top row of frames is for 
frozen potential approximation, the middle shows evolution according 
to the frozen flow approximation and the bottom row is for Zeldovich 
approximation. 

Fig.3 : Same as Fig.2 but for HDM at a = 0.5, 1.0 and 2.0. 

Fig. 7a : An FPA simulation in 2-d. Here we have used a power law 

spectrum with n = — 1. 

Fig. 7b : An N-body frame from a 2-d simulation for the same spec- 
trum as in Fig. 7a. The pancakes are much thinner than the ones in 
the FPA simulation. Also notice the shift of pancakes in this case. 

THESE FIGURES ARE AVAILABLE 
FROM THE AUTHORS ON REQUEST 
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This figure "figl-l.png" is available in "png" format from: 
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